o 

IX. 



o 

CO 



I 

C 

o 
o 



> 
m 
m 

o 
o 



X 



Dynamics of Vesicles in shear and rotational flows: Modal Dynamics and Phase Diagram 
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Despite the recent upsurge of theoretical reduced models for vesicle shape dynamics, comparisons with exper- 
iments have not been accomplished. We review the implications of some of the recently proposed models for 
vesicle dynamics, especially the Tumbling- Trembling domain regions of the phase plane and show that they 
all fail to capture the essential behavior of real vesicles for excess areas, A, greater than 0.4. We emphasize 
new observations of shape harmonics and the role of thermal fluctuations. 
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I. INTRODUCTION 

In the last years there has been an upsurge of inter- 
est in the dynamical response of micro-objects to low 
Reynolds number shear and elongation flows. A large 
number of researchers have dealt with vesicles, includ- 
ing the groups of Steinberg (experiments)^ - — , Gomp- 
per (simulations)^—, Lebedev (theory )2jI2 and Misbah 
(theory and simulations )ii~— . Others have studied, 
along similar lines, their next-of kin, capsules, like 
Seifert, Finken and Kessler (theory and simulations) 15 ' 16 , 
Skotheim and Secomb (theory and simulations )^£ and 
Bagchi and Kalluri (simulations) 18 . A common theme 
is the behavior of the microscopic object at the verge 
of dynamic regime transition. The different investigators 
have discussed an extended phenomenology of intermedi- 
ate and purportedly peculiar motion regimes, supporting 
their findings by the use of the one or the other model, 
or numerical experiment. 

For the simulations, there are two modes of inves- 
tigation: one is direct numerical simulation (DNS) in 
two^ - — or three^ dimensions, which requires compu- 
tational sophistication; the other involves the derivation 
of reduced dynamical models, usually in a perturbative 
framework for nearly- spherical objects, A << 1, which 
yields a number of coupled nonlinear ODEs. The excess 
area of the vesicle A = (A/r^)— 4ir is assumed <C 1, where 
A is the vesicle surface area and ro its characteristic ra- 
dius, obtained from the vesicle volume V = (47r/3)rg). 
Most of the experimental data 1-4 for Tank Treading (TT) 
is in the range 0.05 < A < 2 and for Trembling (TR) 
(called also vacillating-breathing, VB22, and swinging 2 ^) 
is within 0.45 < A < 2. 

The system of ODEs is then studied as a dynamical 
system, either analytically or numerically, and its phase 
space properties are listed. Proposed systems of this sort 
have included two^ r 11 ' 13 , three 2 ^, and recently 14 nonlin- 
ear ODE 12 . A recent review was presented by Vlahovska, 
Podgorski, and Misbah : - which refers to the latest work 
of Misbah and colleagues 1 ^ and includes vesicle results. 
Three important and inter-related issues, often neglected 
by most of the theoretical and numerical work, are: a) the 



effect of thermal fluctuations (because vesicles are small 
and the bending energy of the membrane is comparable 
to the thermal energy); b) the applicability of perturba- 
tive results, to the more readily obtained excess areas, 
A > 0.4; and c) the role of odd modes of the contour 
shape. 

Following the literature, we survey three possible 
regimes of motion, namely (TT), (TR) and tumbling 
(TU). We present here as our benchmark an expanded 
data set including recent experimental worki, with more 
attention to the dynamics of mode interactions. Exper- 
iments of longer duration and higher resolution are on- 
going and will allow us to refine this benchmark. We 
also present a new analysis of our earlier published TT 
datai 2 -, demonstrating the success of the scaling with re- 
spect to the dimensionless parameter A (defined in eq. ((2J) 
below). 

The layout of the paper is as follows: in section [TT] 
we critically review the dynamical models for vesicle 
shapes proposed in literature, which suggest useful non- 
dimensional parameters and provide a basis for laws to 
scale the experimental data. In section IIIII we present 
a new analysis of early and recently published data on 
vesicles and properties of the three dynamical regimes 
observed: in subsection IIII Al we recall the experimen- 
tal setup and the methodology used to analyze vesicle 
contours, and in the following IIII Bt IIII CI and IIII Dl wc 
review experimental findings concerning vesicles in their 
regimes. In section ITVl we consider the broader problem 
of the phase diagram of all regimes of motion. Section M 
presents our conclusions. 



II. DYNAMICAL MODELS FOR VESICLE 
INCLINATION ANGLE AND SHAPE 
EVOLUTION 

A. Models without noise 

Among the reduced dynamical models mentioned in 
the introduction, we refer in detail to three which are 
derived from first principles for a near spherical vesicle, 
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A < 1. A perturbation expansion in powers of v A had 
already been given by Seifert— for Helfrich membrane, 
which is represented by a series of spherical harmonics 
(modes). If harmonics greater than two, thermal noise 
and nonlinear terms above the third order in a free en- 
ergy expansion are neglected, one obtains a system of 
two coupled ODE's: namely the system of Lebedeev, 
Turitsin and Vergeles^ (LTV), that of Danker et al.H 
(DBPVM), and its later variant (KFM). We concen- 
trate on these models because of their central role in the 
recent discussion about regime transitions, because they 
do not contain ad hoc fitted parameters (as the three 
equation system of Noguchi= does), and because they 
can be easily simulated. 

The ODE's describe the vesicle in terms of two dynam- 
ical variables O and ip } associated with the shape and tilt, 
respectively. To avoid confusion we preserve the notation 
of most of the original papers, despite that the tilt -0 will 
be the quantity to be compared to the experimental in- 
clination angle 6, and not O. The different regimes of 
vesicle motion are sought among the attractors of the re- 
sulting dynamical system: for instance a fixed point with 
positive tp is identified with tank treading (TT), a limit 
cycle spanning the whole range [— 7T, tt] for ip is identified 
with tumbling (TU), and a limit cycle spanning just the 
part of the range for ip is identified with trembling (TR) . 
Non-dimensional parameter choices determine the phase 
portrait of the system, the structural stability of phase 
space trajectories, and other issues, which are of concern 
in this paper. Domains in the parameter space leading to 
the one or the other dynamical regime, and their bound- 
ary lines, are referred as to the phase diagram, discussed 
further in section llVl 

The equations of DBPVM read 
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where A = r]i n /r] ou t is the viscosity contrast, s and uj are 
respectively the strain rate and the vorticity of the ambi- 
ent flow, and n is the vesicle membrane bending rigidity 
modulus. We note that while A depends on both A and 
A, Ai and A 2 depend on A alone. In the original formii 
these parameters are expressed using the dimensionless 
shear rate \ (also called the capillary number, Ca), i.e. 
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We remark that this form refers to a pure shear flow, for 
which s = uj = 7/2. The set of parameters {%, A, S, A, A} 
is in fact redundant, at least in pure shear; we include 
here all the definitions for reference, as different sets are 
used in the original papers. 

The DBPVM model differs from the LTV in that it 
introduced "higher-order" terms, with additional param- 
eters Ai and A 2 <C Ai which are multiplied by the small 
parameter \fK in equations ([T]). Further work by the 
same groupie employed a variant of this model (KFM), 
in which only Ai is present. The simpler LTV model, 
in turn, can be obtained from DBPVM just by setting 
Ai = 0, A 2 = 0. 

A further zero-temperature, 14-ODE dynamical model 
was introduced by Farutin, Biben, and Misbah— , in- 
cluding the second and fourth spherical harmonics and 
neglecting all others. This model purports to be most 
realistic, although its results also disagree with the ex- 
periments presented here. The merit of this model seems 
to be its good agreement with unpublished 3D numeri- 
cal simulations. Unfortunately, no details about either 
the analytical or the computational model are given; odd 
modes are absent since they would not be excited on an 
initially ellipsoidal vesicle and the applied shear (both are 
described by j=2 modes), in contradiction to experimen- 
tal observations- -4 . As shown in section UlI CI this is not 
the case, especially in the TU and TR regimes. Also the 
recent direct numerical simulations with thermal noise^ 
mentioned below in section Hi Bl show TR states that are 
greatly distorted and contain odd harmonics. 

The model used by Noguchi and Gompper— which 
also results in a system of two ODE's, can also be quoted 
in this respect, though in part phenomenological. This 
model expresses the dynamics of the vesicle inclination 
angle 9 and of a shape parameter (asphericity) a by 
means of terms which have partially a theoretical jus- 
tification and are partially the result of numerical evalu- 
ations on ellipsoidal shells. The inclination angle 8 and 
the asphericity a of the object are coupled, and ad-hoc fit 
of the free energy function dependent on a is employed. 

To benchmark the different models, aside of results 
available in the literature, we compared the predictions 
about the range of A or A for which the TR regime should 
be observed, at fixed values of A. We choose A or A be- 
cause all proposed models predict TR in a stripe of the 
respective parameter space, with a weaker dependence on 
the second parameter, which is S or \- Table U presents 
a summary of values obtained by different models - some 
as reported by the respective papers, some reproduced 
by our numerical solutions of these models. We con- 
firmed the results for DBPVM, KFM and LTV in a di- 
rect way, integrating the equations in MATLAB using 
the ode45() integrator. We produced phase trajectories 
for selected values of the control parameters, and noted 
the ranges of values of the parameters at which different 
regimes of motion occur. We did not attempt to repro- 
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Paper 


model type 


A 


S 


X 


Atb 


\tr 


Farutin et al^ 


14 ODE, q = 2,4 


0.43 


98.42 


10 


1.83^2.6 


8^ 12 


Vlahovska et al.— 


14 ODE, q = 2,4 


0.43 


12 


1.22 


2 ^ 2.2 


8.91-^9.94 


1.0 


12 


1.22 


2.38^2.9 


8.91-^9.94 


Noguchi and 
Gomppej^^ 


2 ODE, 
phenomenological 


0.44 


12 


1.24 


1.82 -f- 1.95 


7.86 ~ 8.5 


0.91 


36.9 


8 


2.05 ^ 2.2 


5.86-^6.36 


1.44 


23.5 


8 


2.25 ^ 2.4 


4.96-^5.37 


Danker et al.— 
(DBPVM) 


2 ODE, q = 2 
Ai ^ A 2 ^ 
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Kaoui et al.— 
(KFM) 


2 ODE, q = 2, 
Ai / 0, A 2 = 


0.2 


211.61 


10 


1.42 -f- 1.56 


9.33 10.39 


0.43 


98.42 


10 


1.42 ^ 1.57 


5.92-^6.70 


1.0 


42.32 


10 


1.42 ~ 1.615 


3.51-^4.08 


LTV— 


2 ODE, q = 2 
Ai = A 2 = 


0.43 


98.42 


10 


1.41 -j- 1.5 


5.87-^6.33 



Table I. Comparison of ranges of parameters A and A where trembling is observed, according to different models. While we 
relied on information contained in the paper referenced for the first four cases, we numerically simulated the dynamical system 
for the three two-ODE models, DBPVM, KFM and LTV. q denotes the harmonic modes included in the dynamical model. 



duce the results of the 14 equation system of Farutin, 
Biben, and Misbah— (as insufficient details are given), 
and we report here data extracted from their original pa- 
per. The data for the model by Noguchi and Gomppcr 23 
was derived from their Fig. 4b, (equivalent to Fig. 2a of 
Noguchi 26 ). Fig. 2a of Noguchi 26 is in fact a phase di- 
agram for V* = 0.95, 0.9, 0.85 (or A = 0.44, 0.91, 1.44, 
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,2/3 
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which is compared with exper- 
imental data further below in Figure [TSJ We based our 
comparison wherever possible on the parameter values 
A = 0.43 and \ = 10, because these are recurring values 
in the papers we refer to. From Vlahovska, Podgorski, 
and Misbah—, we extracted the range for A from their 
Fig. 6, where S < 12, corresponding to just x = 1-22. In 
all models the boundaries of TR motion are anyway seen 
to depend little on \ a t large S. 

We can draw several conclusions from this compari- 
son. The first and main one is the disagreement in the 
number of independent non-dimensional parameters re- 
quired by each model: all, except LTV, require more 
than two, in contrast with what is observed in the ex- 
periment. The second conclusion is more quantitative: 
none of the models really predicts the observed parame- 
ter range for trembling motion (section IIVI) . If we look 
at A, in particular, the differences in the predictions be- 
tween DBPVM, KFM at different A, and LTV at any 
A turn out to be negligible, but all three models pre- 
dict a narrower range of A for TR than observed. Also, 
the model of Farutin, Biben, and Misbah— gives a pre- 
diction disagreeing with every other model as well with 
the experiment. Finally, we remark that the TR regime is 
inadequately described in these models adopting a frame- 
work of the second and fourth harmonics without thermal 
noise, as a "vacillating-breathing" mode. As was pointed 
in our previous work 3 -^, thermal noise and third harmon- 
ics are crucial for understanding the dynamics of the TR 



state. 



B. Simulations with noise 



Messlinger et al.— present the results of two dimen- 
sional simulations of vesicles based on multi-particle col- 
lision dynamics, or MPC numerical algorithms. This is a 
discrete particle method akin to dissipative particle dy- 
namics and implicitly includes fluctuations due to ther- 
mal motion. It has been developed by Gompper, Noguchi 
and colleagues and applied to simulation of vesicles in 
two and three dimensions. Moreover, Noguchi—, in his 
comprehensive study of vesicles forced with oscillatory 
shear flows, also includes MPC simulations with constant 
shear. His Fig. 1 presents selected shapes during TT, TR, 
and TU. 

It is remarkable that the simulation in 2D, with noise, 
of Messlinger et al4 captures the essence of what we ob- 
serve experimentally. The vesicle shapes, Fig. 3 (and the 
supplementary animation) of their paper show a "swing- 
ing" or "trembling" motion which exhibits concave regions 
of negative curvature in the contour and clearly includes 
higher odd and even harmonics of the radial displace- 
ment. Such shapes are similar to those we show in Fig- 
ures [T] and [T2] below. Also the tumbling motion shown 
by Messlinger et al.— contains higher odd and even har- 
monics, but unfortunately no power spectra which could 
be compared with our experiments ( section 1111 Cp . 

Messlinger et al.— also compare the MPC simulations 
with the two ODE model for the vesicle shape of Noguchi 
and Gompper, mentioned in section Til Al which is stud- 
ied with the addition of noise forcing terms to both equa- 
tions. In Fig. 5, ibidem, the MPC results are compared 
to phase space trajectories of the 2-ODE model, with 
and without forcing (left column MPC, right column, 2- 
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ODE). Trajectories of the a vs. 9 are plotted for various 
parameter values, and it is apparent that, in all cases but 
the first, swinging and tumbling occur intermittently, i.e., 
for a large parameter range the dynamics is strongly de- 
pendent on the noise forcing. This behavior has also been 
remarked in the experiments, and leads to ambiguous 
classification of our TU-TR transition. Unfortunately, a 
complete phase diagram, comparable with others, is not 
included - Fig. 2 of Messlinger et al— has a sample of 
only four points. 

Noguchi and Gompper also discuss in their earlier 
papers 7 * 2 ^, the addition of stochastic forcing terms to their 
system of two ODEs, which are identical in form to those 
in Messlinger et al.—. However, in Noguchi and Gomp- 
per— the effect of forcing is assessed only with respect to 
TT, while all the phase diagrams shown in their other 
papers are obtained in absence of stochastic forcing. 



III. EXPERIMENTAL ANALYSIS OF 
VESICLE DYNAMICS IN SHEAR AND 
GENERAL FLOWS 

In the following, we review experimental findings about 
vesicles in either of the three dynamical regimes men- 
tioned, with particular attention to scaling laws and uni- 
fying parameters suggested by the theories discussed in 
section [TTJ We present a quantitative comparison of the 
old-^ and new^l data on the inclination angles of vesi- 
cles in TT motion, which was obtained by two different 
approaches and analyzed differently in the two sets of pa- 
pers, and comment on the issue of the transition bound- 
ary from TT to TU and TT to TR regimes as function 
of A and A, in old and new sets of the data. Then we 
present a new spectral analysis of shapes, in particular 
for long time series in TR regimes. We also relate to the 
recurrent motion observed in TR and TU, and comment 
about the observed time periodicity. 



A. Experimental techniques, analysis, and 
definition of regime transition lines 

Fluorescent vesicles with prescribed viscosity contrast 
A are prepared and followed when immersed in either 
planar ambient flow with controlled strain s and vorticity 
w, or shear flow with controlled shear rate 7 = 2s in plane 
Couette and channel flow configurations. Many different 
vesicles are imaged in isolation, each for several values 
of either s, u> in plane linear flow, or of 7 in shear flow. 
The experiments have been described previously, and we 
refer to our papers^ — for all details of the procedure. 
We remind that vesicles with 1 < A < 9.19 were studied 
in a shear flow*-— , while only vesicles with A = 1 were 
studied in general linear flow in Deschamps et al.' 4 . The 
bending rigidity modulus of the vesicles was estimated to 
be k — 25ksT in all the experiments. 



Our experimental technique, described in Deschamps, 
Kantsler, and Steinberg—, is capable of determining r 
and A from the full reconstruction of the vesicle in three 
dimensions, but follows the vesicle motion only through 
the imaging of sectional cuts. Care is taken to maintain 
focus and obtain the largest (equatorial) section of the 
vesicle under examination. We maintain that the vari- 
ations in contour shape observed are indicative of the 
dynamics taking part in three dimensions. In Figure Q] 
we show a representative example of an individual vesicle 
(defined by A and A) that changes its regime of motion 
as the ambient general linear flow is changed by variation 

Of Uj/s. 

We analyze the contours of the vesicle images, recon- 
structed at sub-pixel accuracy by an ad hoc fitting al- 
gorithm which locates the brighter edge of the vesicle. 
Such contours are then fitted by the ellipse possessing 
the same tensor of inertia of the vesicle sectional area, 
to define the major axis and the dominant orientation 
— 7r/2 < 9{t) < 7r/2. Two different approaches to deter- 
mine transition thresholds between TT and either TU or 
TR regimes were used in our earlier 3^ and later^ pa- 
pers. In the latter— ii, the regime of motion of the vesicle, 
at time extents during which flow parameters are kept 
constant, is classified automatically according to an em- 
pirical criterion: 



• if M(0(t) - (9)fj > tt/5, then the vesicle regime 
is tumbling (TU); 



• if (9) > M(0(t) - (9)) 2 J, then the vesicle regime 
is tank treading (TT); 

• otherwise, the regime is classified as trembling 
(TR). 

in which (•) denotes the time average. This criterion, 
based solely on the mean and rms fluctuation values of 
the inclination angle and somewhat arbitrary thresholds, 
proved itself simple and robust for isolated vesicles, well 
defined in shape. It works effectively in presence of noisy 
data, and avoids the need of phase-unwrapping the incli- 
nation angle to resolve TU. 

The approach to the data undertaken in the early 
papers 1,2 is somewhat different, and needs to be eluci- 
dated before comparing the various data sets. The early 
experiments were conducted to investigate the TT dy- 
namics, namely the dependence of the inclination angle 
9 on A and A. In order to reduce the scatter in the data 
due to thermal noise, in particular at smaller values of 
A, measurements were averaged over large ensembles of 
more than 500 vesicles of the same value of A (see, for 
example, typical data in Fig. [5]). At that time the value 
of A could only be inferred from the preparation proce- 
dure, and the error on it was found from a representative 
measurement. In the later experiments, instead, vesicles 
were examined individually, and the error on A, amount- 
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Figure 1. Imagery of a single vesicle with X = 1 and A = 1.13 driven into all regimes of motion. Snapshots are recorded every 
1-2 sec. The red line (color online) is the reconstructed sectional contour, further analyzed for harmonics, and the blue dashed 
line is the major axis of the elliptical fit to the contour, according to whose motion the regime is classified automatically. 
The first row shows the vesicle in tank treading (note that the shape is almost elliptical, though random fluctuations are 
appreciable); the second row shows the vesicle during trembling (observe the remarkable changes of shape, with an occasional 
trilobed shape, indicative of a significant third harmonic). The last row shows the vesicle during tumbling. Here the definition 
of 8 is presented. 



ing to about 20%, could be estimated from direct mea- 
surements (see error estimates in Deschamps, Kantsler, 
and Steinberg 3 ). Other parameters, like ro and A, were 
calculated from the cross-section measurements of each 
vesicle with an error of about 20% each, as well. Individ- 
ually, \ depends cubically on ro (eq. ©), and thus can 
vary strongly for different vesicles, even at the same shear 
rate in channel flow. The data ensemble was binned and 
averaged for some class values of A, and was presented as 
(0(A)) for each A available. For example, each curve in 
Fig. 1 of Kantsler and Steinberg 2 includes points result- 
ing from many measurements on vesicles with different 
A and x- At that time, from these ensemble averages, 
we concluded that 9 does not depend on x, in contrast 
to the recent theoretical findings of Kaoui, Farutin, and 
Misbah-i 3 -. 



B. Tank Treading: scaling of 6(A) 

A benchmark for the theories is the prediction of the 
vesicle inclination angle 6 during TT, as function of the 
relevant control parameters. We have addressed the ques- 
tion in our early papers^ 2 -. Very recently Farutin, Biben, 
and Misbah— have analyzed the problem using a part 
of our previously published data (from the figures in 
Kantsler and Steinberg—) and found it in good agree- 
ment with their model. Their presentation can be easily 
compared with that of a larger set of the data, in similar 
variables, shown by Vlahovska and Gracia 27 . Moreover, 
the plot in Farutin, Biben, and Misbah— suggests that 
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Figure 2. Inclination angle from ensemble of single vesicles in 
TT regime as a function of A for A = 1 and r\ — lOmPa s. 
The straight line is is a linear fit to the data points. 



all data refers to x=100, while in Kantsler and Stein- 
berg— the information about x was not provided, for the 
reason just explained in section fill Al To clarify the situ- 
ation we rediscuss the old methodology and present more 
data supporting our view. 

The data for 9 can be plotted versus the one or the 
other representative quantity, like A or A; a functional 
dependence can be sought for single vesicles or for en- 
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Figure 3. Inclination angle from an ensemble of individual 
vesicles in TT regime as a function of A, for two values of A. 
The straight line is a linear fit to the data points. 



semble averages. In Fig. [5] and [3] we show 9(A) for single 
vesicles in the TT regime, respectively for A = 1 and for 
A = 1.8, 2.6. In spite of the scatter of the data due to 
thermal noise, the data collapse remarkably when plot- 
ted against the scaled variable A, and correlate well with 
a linear fit. The value A c , defined by the intercept of 
the fit with 9 — 0, determines the transition, which is 
either to TU or to TR depending on S as discussed fur- 
ther. Compatible A c are obtained in both plots, within 
the experimental error bars. The same was found for old 
data with A = 1 and r\ = l.lmPa s (not shown here). The 
supplementary material at [URL will be inserted by AIP] 
includes tabular data for A = 1.8 and 2.6, with values of 
X for each vesicle. 

In Figure S] we plot the data for (9(A)) of vesicle en- 
sembles. Part of this data were published in Kantsler and 
Steinberg—, and are supplemented here by more values 
of A. All available ensemble averaged data, at different 
A, again collapse when plotted as function of the scaled 
variable A in Fig. [S] (analogously to what was done by 
Vlahovska and Gracia 27 ). The dependence of (9(A)) ap- 
pears as a power law, except for tails at small 9 and 
large A, seen for instance in Fig. 21 in particular for 
A = 3.4, 4.1, 4.9 and 5.3. The tails are strikingly similar 
to those observed in the recent 2D numerical simulations 
of Messlinger et ai— (cfr. their Fig. 4) at comparable 
values of (9) < 0.15 rad and A > 0.7. Their tails are 
explained by the strong amplification of thermal fluctu- 
ations in the vicinity of the transitions to either TR or 
TU. The large scatter of the data at small 9 and large A 
is responsible, in first instance, for deviations from the 
theoretical scaling in Fig. 2 of Vlahovska and Gracia—. 
The extrapolated value A c corresponding to the inter- 
cept 9 — 0, marking the regime transition, cannot in fact 
be compared with a theory that ignores thermal noise. 



Figure 4. Mean inclination angle obtained due to ensemble 
averaging in TT regime as a function of A for vesicles with 
different values A presented on the plot (color online). 



For these tails, the scaling exponent a in the dependence 
A c ~ A a (where A c is defined by extrapolation as with 
A c ), was found in Kantsler and Steinberg— to be about 
— j rather than — i , as predicted later by the theory and 
found for the new data£^. 

Excluding the data points related to the tails at small 
(9), large A and the values of A mentioned above, the 
full set of Fig. [5] can be fitted within the error bars by 
(6>) = 0.81 - 0.46A. This fit provides A c = 1.74 ± 0.2. 
We emphasize that this data was obtained in shear flow 
at x > U which corresponds, according to eq. to 
S > S c = y/3. In light of what became clear later on 
about the phase diagram for vesicles (see section lTvT) . this 
means that the transition at these values of S is from TT 
to TR, and not to TU, as was suggested in Kantsler and 
Steinberg—. 

While there is experimental evidence for scaling of 
(9(A)), a theoretical regression law is more elusive. Theo- 
ries by Seifert— for A = 1, and Misbah.2 2 . and Vlahovska 
and Gracia— for A > 1, provide an exact solution for 
the inclination angle in TT up to its transition to TU, 
9 = \ cos -1 A, giving A c = 1. These theories do not ac- 
count for a transition from TT to TR, whereas the LTV 
theory 9 ^ discusses both possible transitions in detail. 
According to LTV, the transition TT-to-TU occurs for 
S < %/3 at A c = 2/V3 « 1.155, while the transition TT- 
to-TR takes place at A c up to \/2 « 1.41. The solution 
for the inclination angle in TT is found by solving the 
system ([T]) (with Ai = A2 = 0) with l.h.s. equal to zero. 
From the first of eq. ([TJ) , it can immediately be seen that 
this leads to 6 = | cos^ 1 [A cos 0], where can be repre- 
sented in terms of S and A. After some algebra, a closed 
solution 9(A, S) is found, which is weakly dependent on 
S in the range < A < 2/V3, and bounded from above 
by its limit for S = 0, i.e. 9(A) = icos" 1 This 
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solution is displayed, for reference, for the specific values 
S = 10 and 50 in Figs. M and M The LTV solution is 
closer to the experimental data but still disagrees with 
it at A > 0.8 and (9) < 0.35. Moreover the LTV theory 
predicts, for S > V3 and 2/V3 < A < ^2 - 2/S 2 , neg- 
ative vesicle inclination angles, which we didn't observe, 
and a TT motion which is unstable in the third dimen- 
sion and thus in practice not realized. Presenting part of 
the data of Fig. [5] binned in classes of A, rather than A, 
we show clearly in Fig. \§\ that either analytical solution 
matches satisfactorily the data only for (9) > 0.35 and 
A < 1.42. Better agreement between the LTV theory and 
the experiment is found only for the data in the lowest 
A bin, not surprisingly after the basic assumption of the 
theory, A <C 1. Thus, the extension of theoretical results 
beyond A s» 0.15 and small 6, used in particular in the 
recent publications^— , is unreasonable as it violates a 
basic assumption. 

All experimental points corresponding to A=0.15, 
0.24, and 0.42 at all values of (9) (except for two at 
A = 0.15 and (9) close to zero) lie on the fitting straight 
line shown in Fig. [51 within the error bars. As for the 
comparison between the part of the data presented in 
Fig. [5] and the theory of Farutin, Biben, and Misbah—, 
one cannot distinguish the theoretical curves correspond- 
ing to the lower A within the error bars. Only the 
data points shown by the open symbols, corresponding 
to A = 0.77 and 1.42 and 9 < 0.2 and related to the 
enhanced thermal fluctuations at small 9 and large A 
deviate from it. Thus, in our opinion, the quality of 
the experimental data does not allow one to distinguish 
between sets with different A within the error bars. Be- 
sides, as we demonstrated above, the full set of the data 
presented in Fig. [S] is also fitted rather well by a straight 
line in the whole range of (9) (once more, when the data 
points shown by open squares, related to the enhanced 
thermal fluctuations at small 9 and large A at various 
A are excluded). This conclusion brings us back to the 
problem, discussed in the recent theoretical^— and ex- 
perimental papers^, whether a two or three-dimensional 
phase diagram is required to present all vesicle dynamical 
states. According to our statement, only two parameters 
are sufficient to account for the TT data within the error 
bars. If Farutin, Biben, and Misbah 12 were correct, the 
value of the transition A c should depend on A, while we 
claim that for what can be understood from the available 
data, it is not. Further evidence for the regimes of TR 
and TU is discussed below in Section IIV1 

On the other hand, we found surprising and probably 
accidental that the linear approximation for A -C 1 to the 
first solution, i.e. 9 w 7r/4— A/2, describes the data rather 
well inside the error bars, and would give A c = n/2 
1.57. This was already pointed out in our early paper 1 
and actually even used by Vlahovska and Gracia^i. The 
linearization of the solution 0(A, S) of LTV would have a 
slightly milder slope and be as well compatible with the 
data, with A c up to ir/\/3 s» 1.81. 

The new data&i were obtained differently. The de- 
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Figure 5. (9(A)) as a function of A, for the data in TT regime 
presented in Fig. [4] with some typical error bars. The full as 
well as the open squares present all the data with different A 
and A; the open squares indicate the data at small 9 and large 
A at various A, susceptible to enhanced thermal fluctuations. 
The dash-dotted line (dark red online) is the theoretical solu- 
tion 9 = \ cos -1 A; the other three reference lines are the LTV 
fixed point solution discussed in text, respectively for S — 50 
(short dash, blue online), S = 10 (long dash, red online), and 
5 = (dot, purple online). The full straight line is a linear fit 
to the data, based on the full squares only, (6) = 0.81 — 0.46A 
with A c = 1.74. 




Figure 6. (9(A)) as a function of A for the TT data presented 
in Fig. [4] grouped in classes of A. The data points presented 
by open symbols corresponding to A =0.77 and 1.42 and 9 < 
0.2 are related to the enhanced thermal fluctuations at small 
9 and large A. The full straight line is a linear fit to the 
data, based on the full symbols only, (9) — 0.83 — 0.52A with 
A c ~ 1.6 ± 0.2. Reference dashed lines are the same as in 
Fig. [5] (color online). 



7 




0.05 



0.8 1 1.2 1.4 1.6 

A 

Figure 7. {9} versus A in TT regime for vesicles with A = 1 
in four-roll mill experiment. The dashed line is a linear fit 
to the data (8) = 0.83 - 0.56A with A c = 1.49. Data points 
with (6) < 0.1 and standard deviation larger than 0.07 are 
removed. 



termination of the regime of motion of each vesicle as 
a function of A was conducted as explained in subsec- 
tion IIII Al The experiments on vesicle dynamics were 
conducted in two different experimental devices and flow 
configurations. In the plane Couette flow apparatus only 
a minority of vesicles happened to be observed in the 
close vicinity of the transition, since A and A. on which 
the parameter A depends, were not controlled but only 
measured 3 . Besides, a much smaller vesicle population 
was studied than in the old channel flow experiments; 
averaging on 9 was done on short time series for single 
vesicles and not for ensembles. The resulting data is too 
sparse to analyze the TT motion and scaling, and, in 
spite of the smaller error bars and uncertainty in the de- 
termination of A and S, the transition lines on the phase 
diagram in e.g. Fig. 6 of Deschamps, Kantsler, and Stein- 
berg— are marked by rather wide bands. In the four-roll 
mill device^ instead, the control parameter cj/s was var- 
ied for each single vesicle with A = 1 and measured A, 
and more data is available. In Fig. [7]we plot the data for 
(0(t)), averaged over the time versus A obtained in the 
four-roll mill device with the latter procedure, for vesicles 
m TT regime at S > A gain, data can be compared 
to theories only when the noisiest data, i.e. at small 9 
and large A, is excluded. The data that is fairly fitted 
by (9) = 0.83 — 0.56A, giving A c = 1.49, which is compa- 
rable to the value obtained for the old data (Fig. [S] and 
[6]), to the approximate solutions, and even to the upper 
theoretical value A c = y2 provided by LTV for unstable 
TT. 

In summary, the collapse of all experimental data of 
(9 (A)) supports only qualitatively the theoretical sugges- 
tion of the scalin g 9 ' 10 i 22 ' 25 i 27 . This analysis also confirms 
that neglecting the thermal fluctuations leads to the scal- 
ing A c ~ A -1 / 2 , which follows from A c = const, for both 
the old and new experimental data. 




Figure 8. Image of a vesicle in TR, exhibiting multiple out-of- 
focus indentations suggesting a three-dimensional wrinkling- 
like perturbation. Quantitative analysis of the instantaneous 
three-dimensional shape is not possible with our technique. 



C. Trembling: Analysis of experimental vesicle 
contours, dynamics of harmonics and thermal 
noise 

TR, which is the intermediate state between TT and 
TU, turns out to be the key regime to understand the 
vesicle dynamics in a general linear flow. In this regime, 
the inclination angle 9 oscillates around zero. During 
an oscillation cycle a given membrane patch periodically 
experiences both stretching and compression. Because 
of that, the TR dynamics is found to be more complex 
than even TU. The latter is also characterized by the peri- 
odic switching between stretching and compression, but 
the time spent under compression at small inclination 
angles in TR is much longer. This circumstance leads 
to stronger vesicle shape deformations in the TR regime 
due to the volume and surface area constraints and to ex- 
treme sensitivity to thermal noise at small 9. The occur- 
rence of strong shape perturbations and the appearance 
of higher order harmonics resemble very much the wrin- 
kling recently observed and studied in a time-periodic 
elongation flow 28 ' 29 . In the latter case, the control pa- 
rameter, which is the elongation rate, could be varied in 
order to find the onset of the instability and to study the 
nonlinear dynamics of higher order modes, above the on- 
set. Here, like in TR, higher order modes (wrinkles) are 
generated during the compression period, due to the con- 
straints. During compression, strong shape deformations 
present mostly as concavities of the vesicle, producing lo- 
cally a negative surface tension, which in turn initiates 
the instability, with strong sensitivity to thermal noise. 
As the result, both even and odd higher order harmonics 
are generated. Their growth is arrested since compression 
acts for just a brief part of the period, but sometimes, at 
larger noise amplitude, vesicle budding occurs 2 ^. Similar 
effects of vesicle wrinkling and budding and pinching have 
been occasionally observed also in TR (see Figures [FJ [5] 
and [10] and movies [URL will be inserted by AIP] ) . How- 
ever, even the snapshots of a vesicle performing a more 
regular TR (e.g. Fig. [JJ clearly demonstrate a drastic 
difference with the snapshots of the vacillating-breathing 
mode presented in Danker et al.—. In the VB mode 
the vesicle shape remains elliptical, within only even sec- 
ond and fourth harmonics and a vesicle indeed imitates 
breathing. 
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Figure 9. Vesicle with A = 1, A = 2.71 and S = 2.43, generating a protrusion (pinching). The movie is available in the 
Supplementary Material [URL will be inserted by AIP]. 



A = 1.31, s = 0.19/sec, ro = 0.40/sec 

cd o <y cd 
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Figure 10. Trembling vesicle with A = 1, A = 1.72 and S = 9.75, budding. The movie is available in the Supplementary 
Material [URL will be inserted by AIP]. 



To study quantitatively the TR dynamics, we analyzed 
the images of the experiments^, looking at the radial 
amplitude of the contour r(<fi,t), < (j> < 2ir, relative to 
the centroid of the vesicle contour in each image. This 
radial profile is Fourier-decomposed, i.e. r is expressed as 
r((p,t) = ^2 q r q (t)e lq ^ . Observations about the interplay 
of different modes in time are presented in the following. 
We observe as a side remark, that even performing the 
Fourier decomposition in the frame of reference centered 
on the vesicle, there is no a priori geometrical symmetry 
guaranteeing that any particular contour harmonic is null 
or conserved. A translation of the frame of reference 
would alter all modes of the decomposition, but the fact 
that a particular mode (e.g. q = 1) is not null, is not at 
all an indication of a miscentered frame of reference. 

Figure [TT1 displays the time evolution of the inclination 
angle 9(t) and of the amplitude of the lowest harmonic 
modes A q (t) — \r q (t)\ for a typical long trembling se- 
quence, at constant s and u>. Trembling motion is seen 
to be characterized by recurrent, roughly periodic oscilla- 
tions in the amplitude of the zero-th mode, accompanied 
by short duration dips in the amplitude of the second 
mode, correlated with peaks of the third. The latter 
means that the vesicle section periodically departs for 
a short time from its oval shape, to attain a more tri- 
angular appearance. This activity reflects fluctuations 
of the three-dimensional vesicle shape, since the fluctu- 
ations in the observed contour length corresponding to 
the perturbations in Aq can occur only due to 3D effects, 
which are indeed observed in the experiment and resem- 
ble wrinkling (see Fig. [5]). Some instantaneous contours 
of the vesicle during this sequence are shown in Figure 
1121 together with their spectra, notably around the time 
of one of such deformations. The intermediate power 
spectra in Fig. 112b clearly demonstrates the prevalence 



of the third mode over fourth and even second modes at 
the time when the vesicle contour (see Fig. H"2"k ) is sort 
of triangular with an additional concavity. This is ap- 
parent in the movie provided as Supplementary Material 
[URL will be inserted by AIP] . From Figure [TT] we also 
see that the time traces of the first few harmonics are cor- 
related with modes 2 and 3, and that higher harmonics 
are decreasingly smaller and noisier as q increases. The 
inclination angle 9 is seen to oscillate somewhat irregu- 
larly around zero, while the vesicle occasionally performs 
full tumblings. We observe this motion notably in flow 
conditions close to the TR/TU transition (see section ITVl 
below). Such irregularities led us to use the simple and 
robust regime classification criterion mentioned above. 

Figure [T3] provides an example of a vesicle driven into 
the three regimes, evidencing the evolution of the incli- 
nation angle and of the second and third contour modes. 
The second harmonic mode ^(i) peaks and dips irreg- 
ularly during the TR oscillations, with a correlation of 
the dips with the minima of the inclination angle 9; dur- 
ing TU, this mode has minima in correspondence with 
9 ~ (2n + l)n/2, that is twice per full vesicle revolution. 
The increase in As(t) passing from TT to TU and even 
more to TR is apparent, as well as the irregularity of TR 
motion. 

Finally, in Figure [TJJa) we show the spectra of the am- 
plitudes A q = (A q (t)) of the harmonic modes of the con- 
tour averaged in time, for another vesicle with A = 0.64, 
which was also driven in the three different regimes for 
some time. These spectra are typical of all cases ob- 
served, in the following sense: mode 1 is always small; 
modes 2 and 4 have mean amplitudes roughly indepen- 
dent from the regime; mode 3 and 5 are always a few 
times smaller in TT than in TR or TU. The amplitude 
A3 is smaller for vesicles with smaller values of A. These 
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Figure 11. Upper panel: inclination angle 9(t) of of a vesicle 
with A = 1.16 and A = 1, during the course of a long trem- 
bling sequence. Lower panel: amplitudes A q (t) of the first 
harmonic modes of the sectional contour. For convenience of 
plotting, mean values have been subtracted, and time traces 
have been shifted vertically relative one to another 2pm. Ver- 
tical dotted lines are added as a guide to the eye, to mark 
events of occasionahMmWmg (color online). 
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Figure 12. (a) Instantaneous contours in center-of-area coor- 
dinates, and (b) angular spectra of selected snapshots from 
the trembling sequence of figure 1111 Note the significant de- 
crease of A2 and A4 and the corresponding increase of A3 for 
the intermediate contour (color online). 



where a is the surface tension. For large S the last 
term in the denominator can be expressed^ as otq/k = 
0.77SVA. 

For large q, the expression (j4]) can be rewritten as 



rlk B T 



1 9 
ar^q z 



Kq 1 



(5) 



It is remarkable in Figure IT3T b) that the fluctuations for 
TR and TU are quite well described by this latter ex- 
pression, with a — > 0. An explanation of this fact may 
lie in the fact that the effective surface tension, which 
becomes locally negative during the compression of the 
membrane, averages out over an oscillation cycle of the 
vesicle. 



D. Regime transitions and oscillation periods 
of TR and TU 



plots quantify the observation, seen by eye from images 
like those in Figure [IJ that vesicles are "almost ellipti- 
cal" during TT, while they undergo more elaborate shape 
changes during TR and TU. There seems in general to be 
no dependence, either on the regime or on A, of the decay 
spectrum of higher harmonics in the range q ~ 6 -j- 30. 
For the same sequence, the mean squared fluctuations of 
the mode amplitudes, a q = ((A q (t) — (A q )) 2 ), are shown 
in Figure [T4Tb) . It is interesting to note that the fluctu- 
ations of the lower order modes are much smaller in TT 
than in both TR and TU. The spectrum of fluctuations of 
a vesicle can be compared with the theoretical prediction 
of Seifert— in thermal equilibrium 



K{q-l)(q + 2)[q(q + l) + ar%/K 



Figure [T2J as well as the equivalent Fig. 5 of De- 
schamps, Kantsler, and Steinberg— and Fig. 4 of De- 
schamps et al.—, also demonstrate the effect of a quick 
change of the control parameters, 7 in the former and 
lu/s and s in the latter, on the vesicle motion. The time 
series of 6{t) shows an almost instantaneous change of the 
dynamical regime, within about a period, much shorter 
than the characteristic dynamical time r. The theoryii 
and the associated numerical simulations, without ther- 
mal noise, show instead very long transients, of duration 
of the order of x _1 j which are not found in the experi- 
ments. 

The characteristic oscillation period in either TR or 
TU is extracted from all available data like that presented 
in Figures [HI and [T51 To this extent, we determined the 
dominant temporal frequency in A-z(t) using the pburg 
function of Mat lab; this frequency is divided by two in 
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Figure 13. Instantaneous inclination angle 8(t) and amplitudes of the second and third angular harmonics A2(t) and A$(t) for 
a vesicle with A = 0.71 and A = 1 in a general flow, with s and w kept constant for extended time and then quickly changed. 
Vertical (magenta online) dashed lines delimit stretches of constant flow, while horizontal dashed lines mark time averages. 
Note that the transition between one regime of motion and the next is sudden, and the transient appears to be shorter than 
one period (of either TR or TU) (color online). 




Figure 14. (a) Spectra of the mean contour harmonic modes (A q ) for a typical vesicle with A = 0.64 and A = 1, accessing the 
three regimes; (b) spectra of the squared fluctuations 0% for the same sequence. The solid black line is the theoretical form 
e(q, S) given in eq. evaluated for S = 114.11, corresponding to the TT case. The dashed black line is the limiting form 

k k<^ 4 ' ° ' ^ m gh q, the spectra flatten due to the contour reconstruction noise. Horizontal dash-dotted lines corresponding to a 
white noise of amplitude 0.3 and 0.1 pixels are reported for reference, and demonstrate the subpixel accuracy of our analysis 
(color online). 



TU, to account for the observed periodicity within one 
tumbling cycle. The resulting periods T, for the vesicle 
sequences of duration sufficient for a reliable estimate, are 
shown in Fig. 1151 rescaled by either s or uj. Surprisingly, 
we didn't find any correlation between T and r provided 
by eq. ^ . We observed a large dispersion of periods in 
TR, while for large A the periods tend to T ~ 2it/uj in TU 
regime, as expected for rigid body rotation. Correlation 
with A or x, as predicted by theories (e.g. Kaoui, Farutin, 
and Misbah— ) is not observed either. 

Both the the fast regime adjustments and the observed 



oscillation periods in TU and TR disprove, even qualita- 
tively, the picture presented by the reduced theoretical 
models. Instead of a "breathing" of the vesicle shape 
in TR, observed in the models which employ the second 
and eventually the fourth order harmonics and without 
thermal fluctuations, we observe a noisy dynamics. The 
experimental picture found is one of strong mode interac- 
tion and correlation, with a pronounced role of the third 
harmonic, where thermal noise is considerably amplified. 
This underlying mechanism is very different from the 
viewpoint of the papers of Misbah's groupi e 13 ! 22 , though 
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Figure 15. TR and TU periods, rescaled with s and uj (color online). 



the vesicle shape deformations are comparable with those 
presented by Messlinger et al.—. For better appreciation 
the movie of TR dynamics provided in Supplementary 
Material [URL will be inserted by AIP] should be com- 
pared with the snapshots of the VB mode presented in 
Fig. 5 of Danker et al.—. 



IV. PHASE DIAGRAMS - COMPARISON 
WITH RECENT EXPERIMENTAL DATA 

The study of the dynamical systems described in Sec- 
tionllll leads to the construction of phase diagrams. That 
is, regions in the model parameters space, where the same 
regime of motion is attained. A similar approach has 
been applied to capsules^—. While all authors agree 
on the existence of at least the three regimes of motion 
mentioned above, they dissent about the dimension of the 
phase space, and about the position of the regime bound- 
aries in that space. According to LTV, the phase space 
is {S, A} e [0,oo[ x [0, oo [; according to DBPVM and 
KFM the phase space is three dimensional, and better 
represented by the group {Ca, A, A}. The choice of LTV 
scaling has clearly the advantage of simplicity, and, as we 
want to demonstrate, probably accounts for the correct 
scaling in powers of vA, though it agrees only qualita- 
tively with the experiment. To assess this, we plot all 
available data in the one or the other parameter space, 
and compare the results. 

Figure [TBI shows our most recent phase diagram, which 
includes data for vesicles with 0.015 < A < 2.5 from 
the experiments^— and extends the parameter ranges. 
Classification of regimes is automatic for the data of 
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S 

Figure 16. Summary of the motions observed, with coding 
based on the regime - TU, green stars; TR, red triangles, and 
TT, blue circles (full symbols - older datai^; empty symbols 
- newer data^). Data points are plotted with error bars ac- 
counting for the measurement uncertainties. The dash-dotted 
lines are the dividers given by Lebedev, Turitsyn, and Verge- 
les— , while the solid lines are a guide to the eye (for details 
see text) (color online). 



Deschamps et al.— and obtained from the time series 
of the vesicle inclination angle, as explained in Section 
MI CI Data from Deschamps, Kantsler, and Steinberg 3 
had larger error bars and the regime classification was 
done by eye; older data 1 ' 2 is only TT (see supplementary 



12 



material at [URL will be inserted by AIP] for all data in 
tabular form). With this caveat and furthermore, with 
the ambiguities of regime identification of deformed vesi- 
cles close to the transition lines, we claim that cluster- 
ing in different regions is clear. We plot, for reference, 
the dividers of the regimes given by Lebedev, Turitsyn, 
and Vergeles 9 , respectively {5 < V3, A = 2/V3~} for 

TT/TU, 1 5 > \/3, A = for TT/TR, and 

{S > \/3, A = 1.52- 2.12 e- 104S } for TR/TU (best fit 
based on the numerical data of their Figure 9). 

The extent of the experimental regions deviates 
from the model. The solid lines in the figure are 
our eyeball fit to the data, represented by: {0 < 
S < 0.75, A = 1.32} (TT/TU and TT/TR divider) 
and {S > 0.75, A = 1.32 + 0.3 S 35 } (TR/TU divider). 
Note the lower A transition (to tank treading) has zero 
slope on both diagrams in the range of our data, for all A 
sampled. The upper line has finite slope for large S, and 
disagrees with the results of all reduced model simula- 
tions. The divider between TR and TU is also less defined 
on the diagram, in part due to some experimental points 
with large error bars, but mostly because of ambigui- 
ties, namely of vesicles which intermittently flip between 
large-amplitude trembling and full tumbling rotations, as 
exemplified by Fig. [TTJ obtained for S = 4.7 ± 0.3 and 
A = 1.59 ± 0.03. A possible explanation for the fact 
that the upper transition differs greatly from the linear 
theories and the lower transition does not, is that the for- 
mer is associated with a saddle-node bifurcation, while 
the latter with a Hopf type one. It is well known that 
the nonlinear dynamics associated with the saddle-node 
transitions are more sensitive to noise than those near 
Hopf transitions. 

The clustering of regimes would be destroyed if the 
data was plotted, in coordinates {%, A}. In particular, 
all the experimental data for vesicles with A = 1, any 
regime of motion, would collapse on a single horizon- 
tal line. According to Danker et al.—, Kaoui, Farutin, 
and Misbah—, this would happen since the correct de- 
pendence on the third parameter, A, is not taken into 
account. However, contrasting this view we plot a sub- 
set of our data taken in pure shear flow, in Fig. 1171 To 
compare with Fig. 3 of Farutin, Biben, and Misbah 12 , 
we reproduce their transition curves and select among 
the available data vesicles with A close to the values re- 
ported by them. It is rather obvious that the theory^ 2 - 
disagrees with the data even qualitatively. We do not 
attempt comparisons between our data in general flow 
and the Fig. 4 and 5 of Farutin, Biben, and Misbah 12 , 
as not enough data points close to the values displayed 
there are available. 

The dependence on A suggested by Noguchi— (his 
Fig. 2a) for vesicles in uniform shear flow is not matched 
either by the experimental data, as we show in Fig. [18j 
though the discrepancy is smaller than with the theory 
of Farutin, Biben, and Misbah~. The disagreement be- 
tween this model, which includes higher order terms, and 
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Figure 17. Comparison of our experimental results on the 
phase diagram with the theory of Farutin, Biben, and Mis- 
bah— (color online). 



2.5 



< 1.5 



0.5 



.0 




• • ■ ■* • * 1 



. o &c 



o 0.20<A<0.65 ,TT 
v TR 

* 0.70<A<1 .10 ,TT 
▼ TR 

* TU 
■ 1.10<A<1.70 ,TT 

* TR 
TU 

— Noguchi fig. 2a A=0.44 
■ - - Noguchi fig.2a A=0.91 
Noguchi fig. 2a A=1 .44 



0.2 0.5 1 2 3 5 1 2030 50 
S 



Figure 18. Comparison of our experimental results on the 
phase diagram with the model of Noguchi— (color online). 



the phenomenological model of Noguchi—, is also appar- 
ent when comparing Figures [T71 and [T51 

In summary, we conclude that the coordinates of the 
LTV phase diagram (A vs S) look preferable due to its 
simplicity, failing the quantitative agreement with data. 
As we already pointed oufc^, the LTV theory is based 
on some assumptions which are not matched by experi- 
mental observations. However, puzzling as it may be, the 
presentation of a two-parameter phase diagram resulting 
from the self-similar solution of the model provides an 
adequate description of the data. 
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V. CONCLUSIONS 

In attempt to clarify the mechanisms responsible for 
vesicle dynamics, we have carefully reexamined exist- 
ing experimental data, and critically reviewed the re- 
cent theoretical and numerical work. We found that de- 
spite some qualitative features captured by the existing 
reduced models, a good quantitative prediction is not 
achieved. This lack of success may be due to the in- 
complete understanding and modeling of thermal fluctu- 
ations, nonlinear interaction of harmonic modes beyond 
the second (which can form local regions of negative cur- 
vature - which we have seen in the experiments) . Present 
models are derived in the nearly-spherical approximation, 
and even A ~ 0.43 may be considered large. The main 
conclusion of our analysis is that the agreement between 
theory and the experimental observations may not im- 
prove, not even on a qualitative level, by including in the 
existing models just some of the necessary elements. In 
particular, including more even modes without including 
also odd ones and thermal fluctuations, has been proven 
unsuccessful either for the description of the TR dynam- 
ics or the phase diagram presentation. In any event, we 
hope that all future modelers will look critically at the 
effects of noise on transition regimes where (6) < 0.15 
radian and where odd modes of vesicle contours are sig- 
nificant. 
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